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A THEORETICAL AND EXPERIMENTAL INVESTIGATION OF THE LIFT AND DRAG 
CHARACTERISTICS OF HYDROFOILS AT SUBCRITICAL 
AND SUPERCRITICAL SPEEDS ' 

By Kenneth L. W adlin, C’harles L. Shuforu, Jr,, and John R. McGehee 


A theoretical and experimental investigation at subcavitation 
^ speeds was made of the ejffect of the jre e-water surface and rigid 
I boundaries on the lift and drag of an aspect-ratio-W hydrofoil 
I at both subcritical and supercritical speeds and of an aspect- 

I ratio-4 hydrofoil at supercritical speeds. For the aspect-ratio- 

I 10 hydrofoil, tests were made in Langley tank no. 1 and Langley 
[ kink no. 2 at 0.84 <^nd 3.84 chords submergence at subcavitation 

j speeds from 5 to 4^ fps corresponding to Reynolds numbers 

j from 0.18 X 10^ to 1.64 X 10^. For the aspect-ratio-4 

I hydrofoil, tests were made in Langley tank no. 2 at 0.59, 1 .09, 

I 2.09, 3.09, and 4-09 chords submergence at subcavitation speeds 

^ from 15 to 35 fps corresponding to Reynolds numbers from 

I 0.873 X 10^ to 2.04 X 10\ 

( Approximate theoretical solutions for the ejects of the free- 
I water surface and rigid boundaries on lift and drag at super- 

j critical speeds are developed. An approximate theoretical solu- 

1 tion for the effects of these boundaries on drag at subcritical 

I speeds is also presented. The agreement between theory and 

I experiment at both supercritical and subcritical speeds is satis- 

j factory for engineering calculations of hydrofoil characteris- 

tics from aerodynamic data. 

j The experimental investigation indicated no appreciable 
effect of the limiting speed of wave propagation on lift-curve 
slope or angle of zero lift. It also showed that the increase in 
drag as the critical speed is approached from the supercritical 
range is gradual. This result is contrary lo the abrupt in- 
crease at the critical speed predicted by theory. 

INTRODUCTION 

! Airfoils and hydrofoils operate in fluids which differ 
\ principally in density and viscosity, properties that are 

readily treated by the concept of Reynolds number. Since 
such is true, the vast amount of aerodynamic data already 
accumulated becomes available for use in predicting lij^dro- 
foil charactei*istics. The airfoil, however, generally oper- 
ates in an essentially infinite medium, whereas hydrofoil 
applications usually require operation in a limited medium, 
that is, in the proximity of the water surface. Aside from 
the effects of cavitation then, the principal difference be- 
tween airfoil and hydrofoil applications is one of boundaries. 

> Supersedes recently declassified NACA RM L52D23a by Kenneth L. Wadlin, Charles L 
declassified NACA RM L51B13 by Kenneth L. Wadlin, Rudolph E. Fontana, and Charles L 


In restricted areas such as shallow harbors, canals, and 
towing tanks, other boundari(*s are present besides the 
water surface, that is, the bottom and sides. Naturally 
these boundaries also influence the characteristics of a hydro- 
foil, and their effects must be evaluated in order to use aero- 
dynamic data for the prediction of the characteristics of 
hydrofoils under such conditions. 

In addition to the reflective influence of the bottom and 
sides, the finite depth of water limits the speed of propaga- 
tion of the transvei*se waves generated by the hydrofoil. 
This change in flow causes the lift and drag characteristics 
to be different at speeds below this limiting speed or critical 
speed than the}” are above it. 

In the present report available aerodynamic and hydro- 
dynamic theories have been applied to develop an approxi- 
mate method of evaluating the influence of boundaries in 
order to appl}’ existing aerodynamic data to hydrofoils and 
to correct properly data obtained in towing tanks to actual 
open-water conditions. 

Experimental data at subcavitation speeds were obtained 
in two water depths at several depths of submergence at 
subcritical and supercritical speeds and are compared with 
aerodyjiamic data corrected for the boundary effects. The 
boundarj^-correction methods employed are similar to the 
general methods used in wind-tunnel research with the 
additional consideration that the limiting speed of wave 
propagation is taken into account. 

SYMBOLS 

A geometric aspect ratio, 2s/c 

section lift-curve slope at infinite submergence, 
dcijdao 

^2 section lift-curve slope at finite submergence, 

dcifdao 

ai slope of lift curve at infinite submergence, dCLjda 

slope of lift cuive at finite submergence, dCL 2 ld<x 
Co drag coefficient, D/qS 

Coi induced-drag coefficient of rectangular hydrofoil in 

infinite fluid 

^Co^ induced-drag coefficient due to horseshoe-vortex 

images 

Shuford, Jr., and John R. McGehoe, 1952 and contains additional information from recently 
. Shuford, Jr., 1961. 
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induced-drag coefficient due to trailing-vortex 
images 

drag coefficient at infinite submergence, DJqS 
drag coefficient at finite submergence, 
wave-drag coefficient, D^/qS 
lift coefficient, L/qS 

lift coefficient at infinite submergence, LJqS 
lift coefficient at finite submergence, L 2 /qS 
chord of hydrofoil, ft 
section drag coefficient 

section lift coefficient at infinite submergence 
section lift coefficient at finite submergence 
drag, lb 

drag at infinite submergence, lb 
drag at finite submergence, lb 
wave drag, lb 

effective-edge-velocity correction for lift 
Froude number based on depth of hydiofoil sub- 
mergence, V^lgf 

depth of quarter chord of hydrofoil below free- 
water surface, ft 

acceleration due to gravity, ft/sec^ 

depth of water, ft 

lift of hydrofoil, lb 

lift at infinite submergence, lb 

lift at finite submergence, lb 

free-stream dynamic pressure, ^ pV'^ Ib/sq ft 

Keynolds number, Vcjv 
area of hydrofoil, sq ft 
semispan of hydrofoil, ft 
free-stream velocity, ft/sec 

limiting speed of wave propagation or critical 
speed, ft/sec 

induced vertical velocity at three-quarter chord 
due to bound vortex of hydrofoil (surface 
boundary only) 

induced vertical velocity at three-quarter chord 
due to hydrofoil-image bound vortex (surface 
boundar}^ only) 

induced vertical velocity at three-quarter chord 
due to two trailing vortices of hydrofoil (surface 
boundary" only) 

induced vertical velocity at three-quarter chord 
due to two hydrofoil-image trailing vortices 
(surface boundary only) 

induced vertical velocity at three-quarter chord 
due to horseshoe vortex of hydrofoil 
induced vertical velocity at three-quarter chord 
due to hydrofoil-image horseshoe vortex (surface 
boundary only) 

induced vertical velocity at three-quarter chord 
due to hydrofoil-image bound vortices (multiple 
boundaries) 

induced vertical velocity at three-quarter chord 
due to hydrofoil-image trailing vortices (multiple 
boundaries) 

distance of bound vortex measured in free-stream 
direction from three-quarter chord of hydrofoil, 
ft 


y distance to center of image horseshoe vortex, 

measured parallel to lifting line, from center of i 
hydrofoil, ft 

z distance of image bound vortex, measured normal 

to water surface, from li3"drofoil quarter-chord 
point, ft 

a angle of attack, deg j 

«o section angle of attack, deg j 

r circulation strength of vortex, VcC\!2 | 

Ti circulation strength of vortex at infinite submer- j 

gence i 

F2 circulation strength of vortex at finite submergence | 

V kinematic viscosit\^, ft‘/sec 

p mass density, slugs/cu ft * 

a plan-form correction factor for rectangular wings 

(see ref. 7) | 

\l/ h}Tlrofoil submergence parameter (see eq. (16)) i 

DESCRIPTION OF MODELS 

The experimental data were obtained by using 8-inch-chord 
h}^drofoils with aspect ratios of 4 and 10, each supported b\' ^ 

an 8-inch-chord strut intersecting the upper surface of the j 
hydrofoil without fillets. The strut was perpendicular to 
the chord of the hydrofoil. The h^^drofoil and struts were 
made of stainless steel and were polished to a smooth finish 
consistent with wind-tunnel practice. ! 

The hydrofoils had an NACA 64iA412 section which differs | 

from the NACA 64i-412 section 011I3' by elimination of the i 

trailing-edge cusp; the section characteristics of these two | 

are essentially the same (see ref. 1). The strut had an 
NACA 661-012 section. Table I (see page 22) gives the or- 
dinates for the hydrofoil and strut sections as computed 
from references 1 and 2. 1 


APPARATUS AND PROCEDURE 

The tests on the aspect-ratio-10 hydrofoil w^ere made in 
both Langley tank no. 1 and tank no. 2 to obtain two water 
depths. Figure 1 shows a view of the test setup with the 
aspect-ratio- 10 h^'drofoil and the balance attached to the 
structure on the Langley tank no. 2 carriage. The setup 
in Langley tank no. 1 was similar except for the method of 
attachment to the carriage. For the aspect-ratio-4 hydro- 
foil a sting support was used. Figure 2 shows a view of the 
test setup with the hydrofoil supporting sting and the 
balance attached to the structure of the Langley tank no. 2 
carriage. The details of the supporting-sting arrangement 
are shown in figure 3. Figure 4 shows the cross sections of 
the two tanks. Tank no. 1 has a mean depth of 10.64 feet; 
tank no. 2 has a uniform depth of 6.0 feet. 

The hydrofoils were moved vertically by means of a 
motor-driven jacking screw which moved the balance and 
hydrofoil as a unit. Change in angle of attack was obtained 
at the plate attaching the strut to the balance. 

Measurements of lift and drag were made hy means of 
electrical strain gages. The force measurements were made 
at constant speed, angle of attack, and depth of submergence. 
The depth of submergence is defined as the distance from the 
undisturbed water surface to the quarter-chord point on the 
chord line. For the aspect-ratio-10 hydrofoil, tests were 
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I'inuRE 3. — Details of aspect-ratio-4 hydrofoil support arrangement. 
(All dimensions are in inches.) 


Figure 1. — Test setup showing aspect-ratio- 10 hydrofoil and balance 
attached to towing carriage. 
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Tonk no. 2 


Figure 4. — Sectional details of Langley tanks no. I and no. 2. (All 
dimensions are in ft.) 


made at two submergences (0.84 and 3.84 chords) over a 
range of speed from 5 to 45 fps and a range of angle of attack 
from —3.5° to 6.0°. For the aspect-ratio-4 hydrofoil, tests 
were made at five submergences (0.59, 1.09, 2.09, 3.09, and 
4.09 chords) over a range of speed from 15 to 35 fps and a 
range of angle of attack from —3.5° to 4.0°. The change 
in angle of attack due to structural deflection caused by the 
lift and drag forces on the hydrofoil was obtained during the 
calibration of the balance, and the test data were adjusted 
according!}". 

The supporting strut for the aspect-ratio- 10 hydrofoil and 
the supporting sting and strut for the aspect-ratio-4 hydro- 
foil were run alone at the same range of speed, depth, and 
angle of attack as when the hydrofoil was installed. For 
these tests the end of the strut (aspect-ratio- 10 tare tests) 
and the end of the sting (aspect-ratio-4 tare tests) were 
fitted with faired caps. The tares thus obtained were deduct- 
ed from the test data to give the net forces. The net forces 
were converted to the usual aerodynamic lift and drag 
coefficients by using a measured value of p of 1.966 slugs/cu ft 



Figure 2.- Test .^ctnp .'^bowing jisp(.*cl-ratio-4 hydrofoil snsi)en(U*d 
from balance attached to towing carriage. 
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at tlie testing temperatures which were 40° F for the aspect- 
ratio-1 0 tests at 0.84 chord submergence in Langley tank 
no. 1, 44° F for all other tests with aspect-ratio-lO hydrofoil 
in ])oth tanks, and 70° for tests with aspect ratio 4 in tank 
no. 2. All coefRcients were based on the area of the hydro- 
foils. The area of the hydrofoils is 4.44 square feet for the 
aspect-ratio-10 hydrofoil and 1.78 square feet for the aspect- 
ratio-4 hydrofoil. The measured kinematic viscosity of the 
water at the time of the tests in tank no. 1 at 40° F was 
1.85 X 10"^ ftVsec, in tank no. 1 at 44° F was 1.73 X 10”® 
ftVsec, in tank no. 2 at 44° F was 1.83 X 10"^ ftVsec, and in 
tank no. 2 ot 70° F was 1.15 X 10'^ ftVsec. 


EXPERIMENTAL RESULTS 


The basic experimental results corrected for strut deflection 
and drag tares are presented in figure 5 for the aspect-ratio-10 
hydrofoil and in figure 6 for the aspect-ratio-4 hydrofoil as 
curves of lift and drag for each water depth and depth of 
hydrofoil submergence plotted against angle of attack with 
speed as the parameter. The data, converted to coefRcients, 
are presented in figures 7 and 8 in the usual form for aero- 
dynamic data. The strut drag coefficients (based on the 
area of the aspect-ratio-10 h^'drofoil, 4.44 sq ft) plotted 
against speed in figure 9 indicate the range of strut drags 
obtained. 

The lift-curve slopes and angles of zero lift obtained from 
figure 7 for the aspect-ratio-10 hydrofoil are plotted against 
Reynolds number in figure 10. Also included in this figure 
are the corresponding aerodynamic data for the NACA 
653-418 section. These data were taken from reference 3 
and the lift-curve slopes were corrected to aspect ratio 10 
by the equation 


(Xi = 




AE. 


57.3 




( 1 ) 


be masked by Reynolds number effects encountered in the 
tests and by the effects of submergence. 

The variation of drag coefficient with speed for the 10.64- 
foot and the 6.0-foot water depths at lift coefRcients of 0.4 
and 0.6 and depths of submergence of 0.84 and 3.84 chords 
and aerodynamic section drag data at the same lift coefR- 
cients for the NACA 653-418 airfoil section from reference 3 
are shown in figure 11. A comparison of the drag coefficients 
for the two water depths at both lift coefficients and both 
depths of submergence sliows that, with reducing speed, 
when the critical speed in the greater water depth (tank no. 
1, 15.98 chords) was approached, a drag rise occurred 
whereas the drag in the shallower water depth (tank no. 2, 
9 chords) did not rise until its lower critical speed was 
approached. It can be seen that the drag rise increases with 
lift coefficient and decreases with depth of submergence. 
The variation in drag rise with lift coefficient and depth of 
submergence is as predicted by the theory that will be dis- 
cussed subsequently. However, the drag rise was gradual 
rather than the abrupt rise predicted by this theory. The 
trends at the low subcritical spcaals are not too clear since 
they are masked by Reynolds number effects. An indica- 
tion of the possible Reynolds number effects can be obtained 
from the aerodynamic data presented. 

THEORETICAL BOUNDARY CORRECTIONS— SUPERCRITICAL 

GENERAL 

In order to use aerodynamic data for an airfoil in an 
infinite medium to predict t.he cliaracteristics of a hydrofoil 
in the proximity of the water surface and perhaps rigid 
boundaries as would be encountered in shallow water, canals, 
or towing tanks, the influence of these boundaries must be 
evaluated. The boundary condition to be satisfied at the 
free surface is tliat of constant pressure along the surface 
streamlines. The boundary condition to bo satisfied at the 
rigid boundaries is zero normal velocity. 


I 


from reference 4 where Ee is an effcctive-edge-velocity cor- 
rection from reference 5. The hydrofoil data show no 
significant effect of taidv depth at either depth of sub- 
mergence. It is of particular int(*rest to note that, where 
this effect would be expected to be most pronounced, namely, 
in the region between the dashed vertical lines of figure 10 
where the speed in tank no. 1 is subcritical while that in 
tank no. 2 is supercritical, the lift-curve slope and the angle 
of zero lift for a given Reynolds number are essentially the 
same for both tanks. In the region below the critical speeds, 
the trends are not too apparent. The lift-curve slopes de- 
crease and the angles of zero lift increase with decreasing 
Reynolds number, particularly at the shallower depth of 
submergence. Such a tendency is indicated by the corre- 
sponding low Reynolds number aerodynamic data for the 
NACA 653-418 section. The reason for the variation of 
this tendency with depth of submergence is not fully under- 
stood; however, changes in pi*essure distribution due to 
changes in sul)mergence would influence the Reynolds 
niimb(‘r effect. It appears therefore that, if the lift-curve 
sloj)(‘s and angles of zero lift ai'c infliKUiced by the ci'itical 
s|)c'('d, tlu' influ(‘nc(‘ indicatcal by these tests is so small as to 


FREE-SIIRFACE BOUNDARY 

As a first approximation to the three-dimensional problem 
supercritical conditions are assumed, however, with only the 
free-water-surfacc boundary present. The constant-pres- 
sure boundary at tlie free surface can be satisfied by the 
introduction of a horseshoe vortex above the surface which 
has the same direction of rotation as the one which represents 
the loading on the hydrofoil (fig. 12). 

The presence of the image bound vortex does not change 
the direction of the flow relative to the hydrofoil chord line 
in the vicinity of the center of pressure, but it does tend to 
curve the streamlines relative to the hydrofoil chord line. 
The curvature effect is equivalent to introducing camber of 
the hydrofoil in such a manner as to produce a negative lift 
increment. It would seem therefore that a reasonably close 
approximation to the effect of the free surface could be 
obtained by simply evaluating the effect of streamline 
curvature, in addition to the induced-angle effect of the 
trailing vortices, by applying a technique frequently used in 
approximate solutions of aerodynamic problems (see ref. 6). 
This technique involves determination of the circulation F 
required to produce a downward velocity 11^5+ He at the 
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2,800 


X) 2,000 



Angle of ottock, a, deg 


(a) Water depth, 10.64 feet (15.98 chords); depth of submergence, 0.84 chord. 

(b) Water depth, 10.64 feet (15.98 chords); depth of submergence, 3.84 chords. 

Figure 5. — Lift and drag of the aspect-ratio- 10 hydrofoil. 
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3,600 


180 r 



Angle of attack, a, deg 


(c) Water depth, 6,0 feet (9.0 chords) ; depth of submergence, 0.84 chord. 

(d) “U'^ater depth, 6.0 feet (9.0 chords) ; depth of submergence, 3.84 chords. 


Figure 5. — Concluded. 
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Figure 6. — Lift and drag of the aspect-ratio-4 hydrofoil. Water depth, 6.0 feet (9.0 chords). 



• (b) Depth of submergence, 1.09 chords. 
Figure 6. — Continued. 


y63146— 5G 2 


8 


REPORT 1232 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



(c) Depth of submergence, 2.09 chords. 
Figure 6 . — Continued. 



(d) Depth of submergence, 3.09 chords. 
Figure 6. — Continued. 
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(e) Depth of submergence, 4.09 chords. 


Figure 6 . — Concluded. 
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(a) Water depth, 10.64 feet (15.98 chords); depth of submergence, 0.84 chord 
Figure 7. — Lift and drag coefficients of the aspect-ratio-10 hydrofoil. 
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(b) Water depth, 10.64 feet (15.98 chords); depth of submergence, 3.84 chords. 


Figure 7. — Continued. 


three-quarter-chord location which when combined with the 
free-stream velocity V produces a flow tangent to the mean 
camber line of the hydrofoil. Thus, if geometric camber is 
neglected, the hydrofoil angle of attack a is equal to the 
sum of the angles at the three-quarter-chord point induced 
by the hydrofoil vortices and their images located at a 
distance directly above the hydrofoil equal to twice the 
depth of submergence. 

By use of the Biot-Savart law and the notations defined 
in figure 12, the following expressions for the separate con- 
tributions at the line of symmetry to the vertical component 
of the induced velocity at the three-quarter chord were 
obtained: 

The contribution due to the bound vortex of the hydrofoil. 



the contribution due to the image bound vortex, 
r sc 


W,- 


4.[0)‘+(2y)']^(0’+(2/y 


(3) 




the contribution due to the two trailing vortices of the 
hydrofoil. 


^3 = ^ 


r 

2ts 


^)+S^ 


-1 


(4) 


the contribution due to the two image trailing vortices, 


Wr- 


rs 


2x[s*+(2# 




+1 


+s“+(2/)^ 


(5) 


the contribution due to the horseshoe vortex of the hydrofoil, 




( 6 ) 


and the contribution du6 to the image horseshoe vortex 



(7) 
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(c) Water depth, 6.0 feet (9.0 chords); depth of submergence, 0.84 chord. 


Figure 7. — Continued 


By means of equations (6) and (7) a computation of the 
angle of attack a can be made. 

Effect on lift. — In order to estimate the effect of depth of 
submergence on lift, the ratio of the hydrofoil circulation in 
an infinite fluid to that of a hydrofoil at a finite depth of 
submergence for a given angle of attack is obtained; that is, 


r,_ w,+w, , 

(8) 



Therefore (for small angles). 


sc 

r ’ 1^1 

1 ® 



^(2/)^+s^ 


where 




Vca 


and 


r 1 a\ 


Rewriting this equation to get/ in terms of c and substituting 



yields 

di 


1 


A 

1 

+ ^ 1 

, 1 

1+- 


j+<; 


K 

0+^ 


^ (V^^+l + l) 

( 9 ) 


which is the ratio of lift-curve slope at finite depth to that 
at infinite depth when only the free surface is considered. 

For the two-dimensional case, only the induced velocities 
due to the bound vortex W\ and the image bound vortex W 2 
are considered in equation (8) ; thus. 



which is the ratio of the two-dimensional lift-curve slope at 
finite depth to that at infinite depth when only the free 
surface is considered. 

Effect on drag. — In order to estimate the effect of depth 
of submergence on the drag of a finite-span rectangular 
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(d) Water depth, 6.0 feet (9.0 chords) ; depth of submergence, 3.84 chords. 
Figure 7. — Concluded. 


hydrofoil, the drag induced by the hydrofoil images at a 
given angle of attack is obtained from the equation 

( 10 ) 

This relation is not rigorous since it gives an induced drag 
in two-dimensional flow due to the influence of the bound 
vortex at the three-quarter chord. However, for the aspect 
ratios under consideration, when the drag correction is deter- 
mined in the usual manner, that is by evaluating the down- 
wash at the quarter chord, the drag predicted is too low. 
This condition is true even when the spanwise distribution 
of downwash is considered. 

From equation (7) 


where 


Go 

to 

Go 

[l I '' I 




The drag coefficient of a rectangular hydrofoil in an infinite 
fluid is 

The total drag coefficient of a rectangular hydrofoil at a 
given depth of submergence and angle of attack, therefore, is 

( 11 ) 

restricted area 

In order to estim’ate the effect of depth of submergence 
on the lift and drag of a hydrofoil in a restricted area such 
as a shallow harbor, a canal, or a towing tank, a system of 
images (fig. 13) that satisfied the boundary conditions of 
constant pressure at the free-water surface and zero normal 
velocity at the rigid boundaries is required. The boundary- 
induced vertical velocities at the three-quarter chord are 
obtained by computing the combined effect of sufficient 
images to give the desired accuracy. An infinite array of 
images is, of course, required to give an exact value. 
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(a) Lift at all speeds, (b) Drag; depth of submergence, 0.59 chord, 

(c) Drag; depth of submergence, 1.09 chords. (d) Drag; depth of submergence, 2.09 chords, 

(e) Drag; depth of submergence, 3.09 chords. (f) Drag; depth of submergence, 4.09 chords. 

Figure 8 . — Lift and drag coefficients of the aspect-ratio-4 hydrofoil. 
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(a) Tank no. 1. 

(b) Tank no. 2. 

Figure 9. — Variation of strut drag coefficient with speed and depth 
of hydrofoil submergence for the aspect-ratio- 10 hydrofoil. 




Figure 10. — Variation of lift-curve slope and angle of zero lift with 
Reynolds number for the aspect-ratio- 10 hydrofoil. 



(b) Depth of submergence, 3.84 chords. 

Figure 11. — Variation of drag coefficient of aspect-ratio- 10 hydrofoil 
with speed for constant lift coefficient. 



Figure 12. — Horsehoe- vortex arrangement that satisfies the 
free-surface boundary. 


+ 

+ 


LIFT AND DRAG CHARACTERISTICS OF HYDROFOILS AT SUBCRITICAL AND SUPERCRITICAL SPEEDS 


15 
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Figure 13. — Horseshoe- v'ortex arrangement that satisfactorily 
approximates the tank boundary conditions. 


Sufficient accuracy, however, can be obtained with a finite 
array of images. For example, if another row of images 
were added to the top and bottom and another column of 
images to each side of the horseshoe vortex arrangement 
shown in figure 13 (/1=10, tank no. 2, submergence of 0.84 
chord) the additional images would cause a change of less 
than 1 percent in the total induced vertical velocity at the 
three-quarter chord of the hydrofoil. The general equation 
for this velocity for each image vortex (see ref. 7) is: 

For the image bound vortex 

r X r y+g y-s 



and for two image trailing vortices 

W, 


r r y+s r 

' 1 ^ 1 

47r^2“+(y+s)2[_ 

^ ' Vx^+2^+(2/+s)U 


^+{y-sf j- 


where x, y, and 2 define the location of the image with 
respect to the intersection of the quarter-chord line and the 
line of symmetry of the hydrofoil (see fig. 12). 

The ratio a 2 /ai and the drag coefficient Co are obtained 
as previously discussed by substituting W-j+Ws for Wq in 
equations (8) and (10). 

Some results calculated by applying the foregoing theo- 
retical method in tanks no. 1 and no. 2 for estimating the 
effect of submergence on lift-curve slope are shown in figure 
14 for three aspect ratios. 


COMPARISON OF THEORY AND EXPERIMENT 

Lift. — The theoretical results presented in figure 14 are 
compared in figure 15 with the present experimental results 
for hydrofoils of aspect ratios 4 and 10 and with experimental 
results given in references 8 and 9 for hydrofoils of aspect 
ratios 10 and 6, respectively. The ratio a 2 /a-i for the experi- 



(a) A=10; c=0.667 foot. 

(b) A=6; c=0.417 foot. 

(c) .4=4; c=0.667 foot. 


Figure 14. — Effect of depth of submergence on lift-curve slope. 

mental lift-curve slopes for h}^drofoils of aspect ratios 6 and 
10 is the ratio of the lift-curve slope obtained at a given depth 
of submergence to the lift-curve slope (corrected for aspect 
ratio by eq. (1)) as obtained for airfoil data (see refs. 10 and 
11). The ratio ( 12^1 for the experimental lift-curve slopes 
for the aspect-ratio-4 hydrofoil is the ratio of the lift-curve 
slope obtained at a given depth of submergence to the lift- 
curve slope at the greatest depth of submergence. This 
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(a) .4=10; tank no. 1. 

(b) A=10; tank no. 2. (c) ^4=10; tank no. 2 (ref. 8). 

(d) .4=6; tank no. 1 (ref. 9). (e) i4=4; tank no. 2. 

Fioure 15. — Comparison of experimental and theoretical ratio of 
lift-curve slopes for aspect ratios 10, 6, and 4. 

ratio was chosen for the aspect-ratio-4 hydrofoil because the 
experimental lift-curve slope at the greatest depth of sub- 
mergence was approximately 5 percent higher than the lift- 
curve slope (corrected for aspect ratio by eq. (1)) given by 
airfoil data. If the method used for the hydrofoils of aspect 
ratios 6 and 10 had been used, the ratios would be greater 
than 1.0. 

The agreement of the experimental results with results 
given by the theoretical method is generally good. 

Drag. — Results calculated by the restricted-area theoretical 
method for estimating the effect of depth of submergence on 
the drag coefficient are shown in figure 16 for hydrofoils of 
aspect ratios 10, 6, and 4. The magnitude of the increments 
indicates that a correction to airfoil drag coefficients must be 
made to predict Iwdrofoil characteristics at supercritical 
speeds. 

Results calculated by the restricted-area method for both 
tank no. 1 and tank no. 2 are compared in figure 17 with the 
present experimental results for a hydrofoil of aspect ratio 
10. Figures 18, 19, and 20 present similar comparisons for 



Lift coefficient, 

(a) .4=10; tank no. 1. (b) 4=10; tank no. 2. 

(c) 4=6; tank no. 1. (d) 4=4; tank no. 2. 

Figure 16. — Variation of induced-drag coefficient due to the hydrofoil 
vortex images with lift coefficient for hydrofoils of aspect ratios 10, 
6, and 4. 

hydrofoils of aspect ratios 10,. 6, and 4; the experimental data 
for the aspect-ratio- 10 and aspect-ratio-6 hydrofoils were 
obtained from references 8 and 9, respectively, and the ex- 
perimental data for the aspect-ratio-4 hydrofoil are presented 
herein. The agreement of the experimental results with 
results given by the theoretical method is in most cases good. 

THEORETICAL BOUNDARY CORRECTIONS— SUBCRITICAL 

GENERAL 

The speed of propagation of the transverse waves generated 
by the bound vortex of the hyrdofoil is limited to a speed 
which is a function of water depth. This speed is defined 
by where g is the gravitational constant and h is the 
water depth. When the hydrofoil operates below this limit- 
ing or critical speed, the transverse waves travel along with 
the hydrofoil, whereas above this speed the transverse waves 
no longer accompany the hydrofoil. It follows, therefore, 
that the induced effects on lift and drag due to these waves 
are present below critical speeds but not above. The diverg- 
ing waves due to the trailing vortices are not subject to this 
limitation and their effect is present at both subcritical and 
supercritical speeds. The effect, then, of the trailing vortices 
can be computed to a first approximation in the same manner 
at subcritical and supercritical speeds. The effect of the 
bound vortex at subcritical speeds, however, is not the same 
as at supercritical speeds. 
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Tank no. I 


.032 



Tank no. 2 

o Experimental 
Theoretical 

o 



Tank no. I Tonk no. 2 
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\EL L- J ^ 1 \ \ \ I 

0 .2 .4 .6 .8 0 .2 .4 .6 .8 

Lift coefficient, Ci_ 

(a) Depth of submergence, 0.84 chord. 


Figure 17. — Comparison of experimental and theoretical drag 
coefficients for aspect-ratio- 10 hydrofoil. 

Figure 10 indicates that the expected effect of the critical 
speed on lift was either not present or so small as to be masked 
by the Reynolds number effects encountered in the tests and 
by the effects of submergence. It can therefore be assumed 
that to a first approximation the influence of the boundaries 
on lift will be the same as for the supercritical case and that 
only the influence on drag need be considered. 

Since the condition generally encountered in actual appli- 
cations is that of great water depth, most of the theoretical 
work has considered only this case. Mathematical investi- 
gations of the wave drag of a submerged body were made by 
Lamb, who studied the motion of a circular cylinder and a 
spherical body. More exact solutions of these problems 
were given by Havelock, who solved further problems, for 
instance, that of the motion of a submerged ellipsoid. L. N. 
Sretensky (ref. 12) approached the problem of the submerged 
cylinder for both infinite and finite water depths by assuming 
the existence of circulation. Kotchin (ref. 13) gave general 
formulas for the hydrodynamic forces acting on profiles of 
arbitrary shape in water of infinite depth. Vladimirov (ref. 
14) considered the case of a three-dimensional hydrofoil in 
water of infinite depth. Meyer in reference 15 considered a 


J I I L 



Lift coefficient, 

(b) Depth of submergence, 3.84 chords. 

Figure 17. — Concluded. 

two-dimensional hydrofoil in both infinite and finite water 
depths and in reference 16 considered the case of a three- 
dimensional hydrofoil in water of infinite depth. 

DRAG 

In order to estimate the effect of depth of submergence on 
drag, the induced drag due to the hydrofoil trailing-vortex 
images and wave drag must be added to the drag in an 
infinite fluid. 

The boundary-induced-drag coefficient due to the image 
trailing vortices is 

dCn,=C^^(l + <r) (14) 


From equation (5) 

where, for a free surface, 
2s 


V 


Ki- 


V+(2/)* 


2.^sH(2/)^+(|y 
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Lift coefficient, Q 


(a) Depth of submergence, 0.53 chord. 

(b) Depth of submergence, 0.84 chord. 

Cc) Depth of submergence, 1.59 chords. 

(d) Depth of submergence, 2.34 chords. 

(e) Depth of sulirnorgerice, 3.09 chords. 

(f) Depth of subrnergencf!, 3.84 chords. 

Figure 18. — Comparison of expijrimental and thcioretical drag coef- 
ficients for aspect-ratio- 10 hydrofoil (.see ref. 8;. Speed, 25 fps. 


and from equation (13), for a restricted art‘a (tank no. 1 ana 
tank no. 2), 

lu Ft I ^ ~1 

y—s 

2^+(y-sr 


■ 


^/x^+z^+iy 


-S)J 


(15) 


The wave-drag coefficient for a hydrofoil at a given depth 
of submergence and speed is (refs. 13, 15, and 16) 


(16) 

2 — 


where, for a two-dimensional hydrofoil in water of infinite 
depth (refs. 13 and 15) 




\l/=e 





(a) Depth of submergence, 1.0 chord. 

(b) Depth of submergence, 5.0 chords. 

Figure 19. — Comparison of experimental and theoretical drag 
coefficients for aspect-ratio-O hydrofoil (see ref. 9). 


for a two-dimensional hydrofoil in water of finite depth (ref. 
15) 




sinh^[(l-{)eo] 

cosh^C/o-^2 


(18) 


(thep 


arameter Uo is obtained from the relations — r= 


F' tanhUo' 


> 


’ gh Uo 

and for a three-dimensional hydrofoil in water of infinite 
depth (ref. 16) 


where 

/r. ^ \ COS^ ^ COS 20 \ 

(2 COS 0+^cos 6 COS ^ ) j 


0=arc tan 



(17) 
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Lift coefficient, 

(a) Depth of submergence, 0.59 chord. 

(b) Depth of submergence, 1,09 chords. 

(c) Depth of submergence, 2.09 chords. 

(d) Depth of submergence, 3.09 chords. 

(e) Depth of submergence, 4.09 chords. 

Figure 20. — Comparison of experimental and theoretical drag 
coefficients for aspect-ratio-4 hydrofoil. Speed, 25 fps. 


and and are Hankel functions. 

The drag coefficient of a rectangular hydrofoil in an infinite 
fluid is 

C ^ 

CD=Cd+^j^{l + (T) ( 20 ) 


The total drag coefficient of a rectangular hydrofoil at a 
given depth of submergence, angle of attack, and speed is 




l-j-q' I , 1 


ttA 


Stt 


+- 


y2 r 


( 21 ) 


COMPARISON OF THEORY AND EXPERIMENT 

Figures 21 and 22 compare the present experimental results 
for a hydrofoil of aspect ratio 10 with the results calculated 
from equation ( 21 ). The theoretical results were obtained 
by estimating the section drag coefficient q and by adding 


o 

Q 
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(a) Depth of submergence, 0.84 chord. 


Figure 21. — Comparison of experimental and theoretical drag coeffi- 
cients for aspect-ratio- 10 hydrofoil. Water depth, 10.64 feet (tank 
no. 1). 


the calculated induced-drag coefficient of a rectangular hydro- 
foil in an infinite fluid (7^^, the boundary-induced-drag co- 
efficient and the wave-drag coefficient 

The section drag coefficient Cd at low Reynolds number 
was estimated by extending the section drag data of the 
NACA 64i-412 airfoil by comparison (fig. 23) with low 
Reynolds number data for the NACA 653-4 1 8 airfoil sec- 
tion. The boundary-induced-drag coefficient (eq.(14)) 
was obtained by calculating Ki for tank no. 1 and Ki 
for tank no. 2 from equation (15). The wave-drag 
coefficient was computed from equation (16), where the 
values for yp were calculated from equation (17) (infinite 
water depth, two-dimensional hydrofoil), equation (18) 
(finite water depth, two-dimensional hydrofoil), and 
equation (19) (infinite water depth, three-dimensional hydro- 
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(b) Depth of submergence, 3.84 chords. 


Figure 21. — Concluded. 


foil). The values for ^ obtained from equations (17), (18) 
and (19) are compared in figure 24. 

Figures 21 and 22 indicate that the wave-drag coefficient 
for water of infinite depth (two-dimensional hydrofoil) 
added to CD^-\-hCD^ where = gives a better ap- 

proximation of the experimental drag coefficient of a hydro- 
foil at a given depth of submergence and speed than when 
wave drag is calculated for water of finite depth (two- 
dimensional hydrofoil) or for water of infinite depth 
(three-dimensional hydrofoil) . 

This result may be due to the fact that the wave-drag 
theories do not consider both the effect of water depth and 
the three-dimensional case simultaneously, whereas the 
experimental values were at a finite water depth for an aspect- 
ratio-10 hydrofoil. Suitable experimental data for other 
aspect ratios and water depths are not now available to aid 
in clarifying the discrepancy. The difference in the theo- 
retical and experimental results at 5 fps could be an additional 
section drag increment since the section drag coefficient was 
estimated by an arbitrary method. 
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(a) Depth of submergence, 0.84 chord. 


Figure 22. — Comparison of experimental and theoretical drag coeffi- 
cients for aspect-ratio- 10 hydrofoil. Water depth, 6.0 feet (tank 
no. 2). 


CONCLUSIONS 

A comparison of the results calculated by theory and those 
obtained experimentally for hydrofoils having aspect ratios 
of 6 (presented in NACA “WR L-758) and 4 at supercritical 
speeds and an aspect ratio of 10 at subcritical and super- 
critical speeds may be summarized as follows: 

1. A method has been developed which makes it possible 
to calculate at subcavitation speeds, to engineering accuracy, 
the lift and drag characteristics of a hydrofoil from aero- 
dynamic data. The method accounts for the effects of 
submergence of the hydrofoil below the free-water surface, 
the proximity of fixed boundaries, and the limiting speed of 
wave propagation due to limited water depth. 


Section drag coefficient, Drog coefficient, 
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Figure 22. — Concluded. 


Infinite water depth, two-dimensional 
hydrofoil (eq.(l7)) 

Tank no. I, 10.64 ft water depth, two-dimensional 
hydrofoil (eg. (1 8)) 

Tonk no. 2, 6.0 ft water depth, two-dimensional 
hydrofoil (eg. (18)) 



Figure 24, Variation of ^ (see eq. (16)) with speed for infinite water 
depth and for finite water depth at 0.84 and 3.84 chords hydrofoil 
submergence. 


2. There was no appreciable effect of the limiting speed 
of wave propagation on lift-curve slope or angle of zero lift 
at the two depths of submergence investigated. 

3. The increase in drag as the critical speed is approached 
from the supercritical range is gradual. This result is con- 
trary to the abrupt increase at the critical speed predicted 
by theory. 

L.\ngley Aeron.\utical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., Apn7 22 , 1952 . 



Figure 23, — Variation of section drag coefficient with Reynolds 
number (Cl=0.2). 
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TABLE L— ORDINATES OF STRUT AND HYDROFOIL 

(Stations and ordinates are given in inches] 



Upper surface 

Lower surface 

Station 

Ordinate 

Station 

Ordinate 

0. 026 

0.084 

0. 054 

-0.067 

.044 

.104 

.076 

-.079 

.082 

.135 

.118 

-.096 

.179 

.194 

.221 

-.126 

.376 

.279 

.424 

-.164 

.575 

.346 

.625 

-.190 

.775 

.401 

.825 

-.211 

1.176 

.490 

1.224 

-.241 

1.578 

.559 

1.622 

-.261 

1.981 

.611 

2. 019 

-.274 

2.384 

.648 

2.416 

-.281 

2.788 

.673 

2.812 

-.281 

3. 191 

.684 

3.209 

-.275 

3. 595 

.679 

3. 605 

-.259 

3.999 

.661 

4.001 

-.236 

4. 402 

.632 

4.398 

-.207 

4.805 

.592 

4. 795 

-.176 

5.208 

.544 

5.192 

-.142 

5. 610 

.487 

5.590 

-.108 

6.012 

.422 

5.988 

-.076 

6.414 

.349 

6. 386 

-.050 

6.814 

.265 

6. 786 

-.034 

7. 210 

.179 

7.190 

-.022 

7.605 

.090 

7.595 

-.012 

8.000 

.002 

8.000 

-.002 

L.E. radius: 0.083 

Slope of radius through L.E.: 0.168 


Station 

Ordinate 

0 

0 

.040 

.072 

.060 

.087 

.100 

.109 

.200 

.145 

.400 

.200 

.600 

.243 

.800 

.280 

1.200 

.339 

1.600 

.384 

2.000 

.419 

2. 400 

.445 

2.800 

.464 

3.200 

.476 

3.600 

.480 

4.000 

.477 

4.400 

.467 

4.800 

.447 

5.200 

.411 

5.600 

.361 

6.000 

.301 

6.400 

.236 

6.800 

.167 

7.200 

.099 

7.600 

.038 

8.000 

0 

L.E. radius: 0.076 
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